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VLBI imaging of OH absorption: The puzzle of the nuclear 
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ABSTRACT 

Broad hydroxyl (OH) absorption-lines in the 1667 MHz and 1665 MHz transition to- 
wards the central region of NGC 3079 have been observed at high resolution with the 
European VLBI Network (EVN). Velocity fields of two OH absorption components 
were resolved across the unresolved nuclear radio continuum of ~ 10 parsecs. The 
velocity field of the OH absorption close to the systemic velocity shows rotation in 
nearly the same sense as the edge-on galactic-scale molecular disk probed by CO (1-0) 
emission. The velocity field of the blue-shifted OH absorption displays a gradient in al- 
most the opposite direction. The blue-shifted velocity field represents a non-rotational 
component, which may trace an outflow from the nucleus, or material driven and 
shocked by the kiloparsec-scale superbubble. This OH absorption component traces a 
structure that does not support a counter-rotating disk suggested on the basis of the 
neutral hydrogen absorption. 

Key words: galaxies: active - galaxies individual (NGC 3079): ISM - galaxies: 
Seyfert - radio lines: galaxies 



1 INTRODUCTION 

NGC 3079 is an edge-on Sc galaxy housing a LINE R (Low 
Ionisation Narrow Emission Line Region) nucleus jHeckmarJ 
1980) with a long history of observations at various wave- 
lengt hs. The nucleus is also classified as a type 2 Seyfert nu- 
cleus iFord et alJll986l>. The galaxy ha s a systemic velocity 
of 1116 km s" 1 l llrwin fc SeaauisJl99lh . which gives the dis- 
tance of 16 Mpc and hence 1 milliarcsec (mas) corresponds 
to approximately 0.08 pc, adopting Ho = 75 km s _1 Mpc - . 
Observations of optical emissions with the Hubble Space 
Telescope (HST) have shown a number of outflowing fila- 
ments in a 'superbubble', produced by a combination of stel- 
lar winds and supernova explosions in a site of massive star 
formation. X-ray observations with Chandra show a clear 
correspondence of Ha-line filaments with thos e in X-rays 
at th e distance of ~ 65 pc from the nucleus JCecil et alJ 
2002). Throughout several X-ray observations, a highly ob- 
scured active galactic nucleus ( AGN) and iron lin e emission 
towards the nucleus were found JCecil et al]l2002fr . The cen- 
ter of the galaxy has been resolved by radio interferometry 
and reveals a double radio lobe extending 1.5 kpc on either 
side of the major axis of t he galaxy and a puzzlin g nuclear 
structure on parsec scales lllrwin fc Sea quist 1988J). Due to 
the complexity of the radio structure, the interpretation of 
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the location of a true nucleus is not straightforward. 
Broad (200-400 km s _1 ) and deep (r ~ 0.5) neutral hy- 
drogen HI absorption and weaker but equally broad and 
deep OH absorption are found towards those compact radio 
sources in the nucleus. These absorptions have been studied 
in detail at various angular resolutions by radio interferom- 
eters. Based on the Very Large Array (VL A) observations 
at 1 arcsec resolution, iBaan fc IrwirJ l|l995|) considered that 
the two components of HI and OH absorptions are associ- 
ated with an obscuring torus in the nuclear region. Multi- 
Element Radio Linked Interferometer Network (MERLIN) 
observations at sub-arcseco nd resolution resol ved the rota- 
tion of the HI absorption JPedlar et alj fl996') in the same 
sense as the rotational trend trace d by HI emission and 
CO(l-O) emission o n galactic scales lllrwin fc SeaauisJi99lt 
llrwin fc Sofuell 992'l . According to Very Long Baseline Inter- 
ferometry (VLBI) observations at an angular resolution of 
10-15 mas, three HI absorption components viewed against 
the resolved double-peaked radio source appear to be ev- 
idence for co unter rotation relative to the rotation in the 
outer galaxy ( Sawada-S atoh et al ] |200dh 
The strong and highly intensity-variabl e water maser emis- 
sion is known to exist in NGC 3079 (iHenkel et alJ Il984t 
lHaschick fc Baarj|l985T) . Earlier attempts for resolving the 
maser emis sion have been mad e since the late 1980s us- 
ing VLBI (lHaschick et al]|l99Pl) . At present, it is under- 
stood that the H 2 maser components are distributed 
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Table 1. Properties of NGC 3079. Coordinates in this table were 
phase centre values used for the data correlation at JIVE. Veloc- 
ity with respect to LSR for radio definition is adopted from Irwin 
& Seaquist (1991). Optical classification was made by Heckman 
(1980) and Ford et al. (1986). Radio flux densities were measured 
by the VLBA (Sawada-Satoh et al. 2000). Infra-red (IR) lumi- 
nosity is from Baan (1989). X-ray luminosity is from lCecil et all 
i2002h . We assume that H = 75 km s _1 Mpc _1 . 



RA (J2000) 


10 h 01 m 57.8050 s 


Decl (J2000) 


55°40'47.080" 


Systemic Velocity (21-cm HI) 


1116 ± 1 km s" 1 


Distance 


16 Mpc 


Inclination (optical) 


84° 


Optical class 


LINER/Seyfert 2 


F„ (1.4 GHz) 


5.7 ± 1.7 mjy 


(8.4 GHz) 


14.2 ± 0.7 mjy 


(15 GHz) 


24.4 ± 0.6 mjy 


(22 GHz) 


12.1 ± 0.7 mjy 


Lir 


3 10 10 erg s" 1 


Lx-ray(0.1-6.5keV) 


1.0 10 42 erg s" 1 


Lx-ray(Pc linc;6.4kcV) 


3.3 10 41 erg s" 1 



within an elongated (dis k-like) structure wi thin a few par- 
sees from the nucleus jTrotter et all Il998t) . Most of the 
maser components detected with VLBI lie in a velocity 
range Vlsr (LSR refers Local Standard of Rest.) = 956- 
1190 km s _1 and are not asso ciated with any of the .je t 
components llTrotter et al.ll99<l : ISawada-Satoh et al.l200C() . 
They are distributed nearly north-to-south, and are aligned 
roughly with the molecular disk traced by CO (1-0) (Position 
Angl e (P. A.) = 15°) with a radius of 750 pc (llrwin fc Sofud 
1992). Several new weak components were recently detected, 
which e xtend the velocity ran ge to about Vlsr = 1350 
km s _1 IIHaeiwara et alJl2002al) . Given the velocity distri- 
bution of all the maser features observed to date, the maser 
spectrum shows a rough symmetry w.r.t th e systemic veloc- 
ity of the galaxy (Vlsr — 1116 km s -1 ) ( Haei wara et alJ 
20023). It has been argued that there is a rotating parsec- 
scale molecular gas disk or torus obscuring an AGN in 
NGC 3079, but no concrete evidence has been presented in 
the literature like the one in NGC 4258 jHerrnstein et alJ 
119981) . A recent study bv lMiddelberg et all (120031) presents 
evidence for five distinct radio continuum sources in the nu- 
cleus of NGC 3079, complicating the interpretation of its 
puzzling structure. Besides the three components A, B, and 
C that are detected at higher frequencies and are seen at 22 
GHz together with the H2O masers, there are two other com- 
ponents E and F that are visible only at lower freq u encies , 
following the convention used in llrwin fc SeacmisH (l988). 
These additional components have the same orientation and 
are equally spaced as the 22 GHz components A and B, but 
they are shifted about 25 mas ~ 2 pc to the east; they form 
the 18-cm radio signature that is observed in HI and OH 
absorption experiments. 

In this paper, we present VLBI observations of OH ab- 
sorption at the nuclear region of NGC 3079. With observa- 
tions using European VLBI Network (EVN) at the highest 
spectral-line sensitivity of any existing VLBI facility, we aim 
to review the kinematics of the circumnuclear region in the 
galaxy. 



2 OBSERVATIONS 

The central region of NGC 3079 was observed at 18-cm on 
2000 November 15 with the EVN, which consisted of eight 
telescopes: Cambridge, Effelsberg, Jodrell Bank (Lovell), 
Medicina, Noto, Onsala, Torun, and the Westerbork 
phased-array. We observed the 1667 and 1665 MHz main 
line transitions of the 2 f7 3/ / 2 , J=3/2 ground-state of OH 
in absorption toward the central radio continuum source in 
NGC 3079. 

NGC 3079 was observed in a phase-referencing mode over 
a period of 5.7 hrs, interspersed with observations of a 
phase calibrator source J 0957+5522, which is located 
about 1 degree from the target source. The total time 
of one observing cycle was 13 minutes with 10 minutes 
for NGC 3079 and 3 minutes for J 0957+5522. The total 
time spent on the target was about 3 hours. The data 
were recorded in both left and right circular polarizations 
using a single intermediate frequency (IF) band with 8 
MHz bandwidth. The IF baseband was subdivided into 128 
spectral points, yielding a frequency resolution of 62.5 kHz, 
or 11.2 km s in velocity at the source distance. The IF 
velocity coverage is 1440.5 km s~ x . In order to cover both 
the main line transitions, we centered the IF on Vlsr = 
1350 km s _1 and the Doppler velocity center was referenced 
to the 1667 MHz line. Hereafter, all velocities are in the 
radio convention and with respect to LSR. The system 
temperature and antenna sensitivity for each EVN antenna 
over the observing frequency range from 30 K — 100 K dur- 
ing the observing run. DA 193 was measured in the middle 
of the observations for calibrating the absolute amplitude 
gain and the bandpass correction. The data-recording rate 
was 128 Mbit per second with a MklV terminal. 

The correlation of the data was performed at the EVN 
MklV correlator at the Joint Institute for VLBI in Eu- 
rope (JIVE). Data analysis was made using the NRAO 
AIPS package. After the delay and delay rate calibra- 
tion using J 0957+5522, the bandpass corrections were ap- 
plied. We discarded the data from six baselines, Noto(NT)- 
Medicina(Mc), Mc-Torun(Tr), Tr-NT, Cambridge(Cm)-Tr, 
Cm-NT, and Cm-Mc baseline, due to inadequate data qual- 
ity. The continuum visibility data set was generated from 
the spectral-line visibility data set by averaging absorption- 
free channels using the AIPS task of UVLIN. The contin- 
uum visibility data set was then used for self-calibration to 
improve the image sensitivity. After iterations of the self- 
calibration in both phase and amplitude, the solutions were 
transferred to the spectral-line visibility data. The contin- 
uum and spectral-line visibilities were imaged employing 
IMAGR. 



3 RESULTS 

Fig. 1 (left) shows the EVN spectrum of OH absorption 
integrated over the nuclear radio continuum source (Fig. 4) 
with natural weighting with a spectral resolution of 11.2 
km s _1 . The systemic velocity of the galaxy Vlsr = 1116 
km s _1 and the peak velocities of HI absorption and water 
maser emission are denoted. The rms noise level is 0.7 mjy 
beam -1 per channel. The OH main-line transitions of 1667 
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MHz and 1665 MHz are clearly detected against an unre- 
solved continuum source at the centre of NGC 3079. Both 
transitions have two distinct components, one of which is 
blue-shifted (OH1) and the other is nearly centred on the 
systemic velocity (O H2). This is consistent with the spectra 
at lower resolutions l|Baan fc Irwinl Il995l) . By comparison 
with the single-dish spectrum, the OH absorption lines ob- 
served in these observations have lost their broad wings. The 
integrated intensities of the 1667 MHz OH absorption esti- 
mated from the single-dish and t he EVN spectrum in Fig. 
1 (left) are about 8.3 Jy km s -1 (Ba an fc Irwinlll995l) and 
0.74 Jy km s - , respectively. Only 8.9 percent of the ab- 
sorption intensity was recovered in the EVN observations. 
This accounts for the narrower line-widths of the OH ab- 
sorption spectra obtained with the EVN (Table 2). The OH 
emissions in the wing of the absorption at V lsr = 910 and 
aroun d Vlsr = 1260 km s -1 suggested bv iBaan fc Irwiiil 
l|l995h coincide with weak emission features seen in our EVN 
spectrum, although both of them need to be confirmed in 
further observations. Two Gaussian components have been 
fitted to the absorption profiles of each of the transition (Fig. 

1 right); the results are listed in Table 2. The resultant 
opacity and column density are listed in Table 3 together 
with those of the previous VLA and VLBI observations ob- 
tained at different angular scales. The OH absorption veloc- 
ity center (Gaussian-fitted) is Vlsr (OHl)= 1011.9 ± 0.9 
km s -1 and Vlsr (OH2)= 1113.5 ± 2.0 km s -1 , while the 
centre velocities derived from the VL A-A data are Vlsr = 
1013 km s" 1 and Vlsr = 1114 km s" 1 iBaan fc Irwiii l995l. 
These values are consistent within the spectral resolution 
of one channel. The HI absorptions are pea ked at Vlsr = 
1010 km s" 1 , 1120 km s" 1 , and 1265 km s" 1 iBaan fc Irwinl 
1995), where the red-shifted third component has no OH 
counterpart. The ratios of the double peaks in both transi- 
tions are quite consistent, which may imply that both main- 
line features are associated with the same continuum source. 
The line ratio 1667 to 1665 MHz in the Gaussian-fitted pro- 
files is 1.5, which is lower than the Local Thermodynamic 
Equilibrium (LTE) value of 1.8. 

Velocity-integrated maps of the 1667 MHz OH absorption 
integrated over 11.2 km s -1 intervals and the CLEAN 
map of a 18-cm radio continuum source are shown in Fig. 

2 (All are uniformly weighted.). Fifteen velocity channel 
maps cover the velocity range of the absorption from 990.4 
km s -1 to 1147.7 km s -1 . 

The continuum emission was not resolved in this experiment. 
The peak flux density of the continuum map produced with 
uniform weighting is 8.8 mjy beam -1 and the rms noise level 
is 0.25 mjy beam -1 . The peak flux density and the noise 
level of the continuum map produced with natural weighting 
are 11.8 mjy beam -1 and 0.085 mjy beam -1 , respectively. 
The integrated flux density of the continuum is 14.3 ± 0.7 
mjy, while the 21-cm integrated flux densi ty obtained b y 
single-dish measurements is 760 ±31 mjy (ICondonl ll983). 
Hence, more than 95 percent of the continuum flux is miss- 
ing in this EVN experiment. Because the limited angular 
resolution of our EVN observations compared with the past 
VLBI observations at higher frequencies, the nuclear con- 
tinuum remains unresolved without showing any structure. 
The EVN synthesized beam is larger than the angular sep- 
aration between the various nuclear components (Fig. 7), 
which limits the discussion on the possible associations of 



the OH absorbing gas with individual sources in later sec- 
tions. 

A position-velocity (PV) diagram of the 1667 MHz absorp- 
tion is presented in Fig. 3. A cut was made through the 
continuum emission at P.A. = 15° (rotated clockwise by 15 
degrees), which aligns a c ut along a major ax i s of the ro- 
tatio n traced by CO(1-0) iSofue fc IrwirJll992l : Ikoda et alJ 
20021) and which was also used from the PV diagrams in 
Baan fc Irwinl (Il995t) . The gradients in the 1665 MHz line 
show the same trend as those seen in the 1667 MHz line. 
Note that there are significant position offsets of ~7-8 mas 
in declination between the two peaks of OH1 (blue-shifted) 
and OH2 (syste mic), which could n ot be measured in the 
PV diagrams of IBaan fc Irwinl il995l) at 1 arcsec resolution 
with the VLA-A. The velocity contours of OH1 and OH2 
ar e separate, as c ompa red with the PV diagrams in Fig. 7 
in lBaan fc Irwinl dl995l) . 

Fig. 4 shows a naturally weighted 18-cm continuum map 
superposed on the integrated 1667 MHz OH intensity (0th 
moment) map. The absorption is apparently concentrated 
towards the centre with a weak outward extension. This is 
similar to the east- west elongation of the HI and OH absorb- 
ing gas seen at 1 arcsec resolution and caused by the orien- 
tation and separation (~ 25 mas or ~ 2 pc) of the L-band 
components E and F <|l rwin fc Seaauisdll988l:lBaan fc Irwinl 
ll99alMiddelberg et alJl2003l) . 

Fig. 5 displays the OH velocity field of the two components 
(1st moment). The velocity field of OH1 varies from south- 
west to northwest in P.A.= 60°, while that of OH2 is seen 
approximately from north to south in P.A.= 145°. The ve- 
locity gradients of OH1 and OH2 are ~ 10 km s -1 pc -1 
and ~ 13 km s -1 pc -1 , respectively. The kpc-scale CO(l- 
0) velocity gradient is found to be also in the north-to- 
south directio n and has a value of 0.85 km s -1 pc -1 and 
a P.A.= 15° iSofue fc Irwinl Il992l) . The compact nuclear 
CO (1-0) core at PA= 4°, tilted by about 10° w.r.t the 
kpc-scale CO(1-0) disk, within central 125 pc s hows a rigid 
rotat i on with a velocit y of ~ 300 km s -1 iSofue et alJ 
l200lt iKoda et alJ 1200211 . This yields the veloc i ty gra dient 
of 2.4 km s -1 pc -1 . Jus t as IBaan fc IrwirJ (jl995l) also 
ISawada-Satoh et alJ feOOd) have interpreted the equivalent 
HI components to originate against the two continuum com- 
ponents A and B, which in reality were components E and F 
ilMiddelberg et alJl2003ft . The velocity differences estimated 
from the CO velocity gradients between the two contin- 
uum components F (A) and E (B), separated by 25 mas 
or 2 pc, are only 1.7 km s -1 or 4.8 km s -1 . Such values 
are not consistent with the observed velocity difference of 
101.6 km s -1 between OH1 and OH2, which would suggest 
that the two detected OH absorptions arise from neither the 
kpc-scale CO disk nor the nuclear CO core. MERLIN ob- 
servations at ~ 400 mas resoluti on resolved the ro tational 
motion of the HI absorbing gas dPedlar et alJll996ft . which 
matches the sense of the rotation traced by the HI emission 
Jlrwin fc Seaauisdll99lh . The directions of the HI velocity 
fields in the three HI VLBI absorption c o mpon ents have not 
been presented in ISawada-Satoh et alJ <l2000l) . which pre- 
vents a comparison with the OH velocity fields in Fig. 5. 
Comparisons of the intensity distribution of OH1 and OH2 
for each transition are presented in Fig. 6. The intensity 
maps were produced from the uniformly weighted spectral- 
line cubes. We find that there is a difference in the spatial 
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distribution in both 1667 MHz and 1665 MHz absorption, 
which is even more distinct in the 1665 MHz line. The dis- 
tribution of the blue-shifted 0H1 is downward w.r.t that of 
the systemic OH2 by approximately 8 mas, or 0.6 pc in dec- 
lination. This trend is also true in the PV diagram of Fig. 
3, where the difference is clearly seen for both main-lines. 



4 DISCUSSION 

4.1 The puzzle of the nuclear radio continuum 

NGC 3079 is known to exhibit a complex radio structure 
in the nuclear region. Recent VLBI studies reveal a mul- 
tiple milli-arcsec-scale nuclear so urce structure with fiv e 
frequency dependent components jMiddelberg et al\ \2003). 
The continuum components A, B and C are only present 
in maps at 22 GHz but are totally obscured at L-band. 
The identification of a true nucleus in the galaxy has been 
a subject to debate. If the H2O maser structure is indeed 
confirmed to be a disk with a rotation axis misaligned 
by 15-20° with respect to the kpc-scale CO(l-O) disk 
at P.A.= 30°, then the nucleus would lie close to the 
A-C-B line (Fig. 7). iTrotter et all <1998ft proposed that 
the dynamical center lies between components A and B, 
more specifically at the intersection of the jet axis and a 
north-south distribution of the H2O masers. The north- 
south distrib ution of the H2O m aser was more convincingly 
confirmed in iKondratkol J2003T) by measuring the positions 
of the new red-shifted H2O maser features reported in 
lHagiwara et alJ i2002al) . An H2O disk centered close to C 
and having a systemic velocity of 1120 km s _1 would be 
consistent with the H2O components at 1123 km s _1 at 
C, the systemic velocity of the HI absorption of 1116 
km s _1 , and the systemic OH component (OH2) at 1114 
km s _1 . The models based on the velocity sign ature of 
the OH and HI absor ption feaan fc Irwinl Il99^ and of 
the HI VLBI signature l|Sawada-Satoh et al l2000t) refer to 
components E and F and do not reflect the location of the 
nucleus. Various monitoring programs have been executed 
to detect a velocity drift of the systemic components as 
seen in NGC 4258, but they have failed to find any drifting 
component to confirm the disk signature (e.g. Baan & 
Haschick 1996; Hagiwara et al. 2002a). 

Considering the available evidence, the true nucleus is 
hidden at L-band and lies close to C on the line connecting 
components C and B. The nature of the L-band components 
E and F is not yet certain. These components could be 
circum-nuclear starburst regions but for that E has too 
high a brightness temperature and could be a supernova 
remnant. 



4.2 The nature of the OH absorption 

An important result of our EVN observations is the inter- 
pretation of the two different velocity fields traced by the 
two OH absorption components on scales of 10 mas towards 
the unresolved radio continuum nucleus. The spectrum in 
Fig. 1 shows that most of the OH absorption has been re- 
solved out in our data set, as compared with single-dish and 
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Figure 1. left: Hydroxyl absorption spectrum of the 1667 MHz 
and 1665 MHz transitions, obtained with the EVN in November 
2000. An rms noise per channel is 0.7 mjy beam -1 . The velocity 
resolution is 11.2 km s . Velocities in the spectra are scaled 
in the radio LSR convention. The adopted systemic velocity of 
NGC3079 is V L SR = 1H6 km s -1 . The peak velocities of HI 
absorption and H2O maser are denoted by arrows. The results 
of 2-D Gaussian- fitting of the spectrum are listed in Table 1. 
right: Gaussian models fitted to the OH spectrum of the 1667 
and 1665 MHz transition. The velocity is referenced to Vlsr = 
1013 km s , the Gaussian-fitted peak velocity of the 1667 MHz 
OH1 in Table 2. 
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Figure 2. Velocity channel maps of uniformly weighted 
continuum-subtracted spectral line cube where the OH absorption 
at 1667 MHz is present. The contour levels are -8,-7,-6,-5,-4,-3, 3 X 
0.8 mjy beam - x (l a). 18-cm uniformly weighted continuum map 
is also shown with the contour levels at -5,5,10,15,20,25,30,35,40 
X 0.25 mjy beam - 1 (1 a). The center (0,0) position is 
R.A.(J2000)=10 h 01 m 57 s .805 Dec.(J2000)=+55°40'47".08. The 
synthesized beams are plotted at the top-left and the bottom- 
right corner. 
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Figure 4. 18-cm naturally weighted continuum map (contour) 
superposed with a 1667 MHz OH intensity (0th moment) map 
(grey-scaled in Jy beam - m s — 1 ). The peak flux density of 
the continuum is 10.3 mjy beam . The contours are plot- 
ted at -5,5,10,20,30,40,50,60,70,80,90,100,110,120 x 0.085 mjy 
beam - x (l cr). The synthesized beam (70 mas x 39 mas, P.A.= 
—72°) is plotted on the lower left corner. 
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Figure 3. Position-velocity (PV) diagram produced from the nat- 
urally weighted spectral-line visibility cube, sliced along P.A.= 15 
degre es to align the major axis of the CO rotation JSofue &: Irwinl 
1992) and the cuts used in PV diagrams in (Baan & Irwin 1995). 
The contours are plotted at 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 
percent of the peak intensity of -6 mjy beam -1 . The positive 
contour corresponds to +0.6 mjy beam -1 . 



VLA spectra in lBaan fc Irwinl . This implies that we 

observe compact gas that extends no more than the beam 
size of 45 mas, or ~ 4 pc, which is almost equal to the size 
of the background continuum. The OH absorbing gas has 
been interpreted primarily to be a p art of circum-nuclea r 
gas traced by the CO(1-0) emission feaan fc Irwinlll995l) . 
A single black body fit based on the measurement of far- 
infrared (FIR) flux densities at 12 - 100 urn in NGC 3079 
gives a temperature of 43 K and the maximum extent of 
this large-scale FIR source (e.g. torus) of about 130 pc. Our 
EVN data do not provide compelling evidence for such a 
large-scale torus. 

4-2.1 The absorbing gas structure 

The nucleus of NGC 3079 hosts a LINER or a type 2 Seyfert 
nucleus. Therefore, the nucleus could be obscured by an 
edge-on dusty torus or intervening medium along the line 
of sight. Because of the systematic sim ilarity of the HI and 
OH absorptions, iBaan fc Irwinl (Il995t) suggested a possible 
association of the two HI and OH absorbing gas components 
with the double continuum source (Fig. 7). With a 1.0 
arcsec (corresponding to 80 pc ) resolution of the V LA in 
A Configuration, the data in IBaan fc Irwinl il995tl were 
insufficient to extract precise connections between these 
absorption components and the nuclear radio continuum 
structure. 

Interp retation of the VLBI data in ISawada-Satoh et alJ 
(2000) shows that the three HI components are resolved 
and that the column densities of each component against F 
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Figure 5. Mean velocity field (1st moment) of the two 1667 
MHz OH absorption components. Contours are plotted every 10 
km s" 1 in LSR velocities with the gray scales indicated on the 
top. left: blue-shifted component (OH1), right: systemic velocity 
component (OH2). 



(blue-shifted HI), F+E (systemic HI), and E (red-shifted 
HI) are almost the same, suggesting that E and F are 
uniformly obscured by the ISM foreground to the nuclear 
sources. Our EVN data failed to resolve E from F, so that 
we cannot identify the individual contributions for OH ab- 
sorption and th e OH column densitie s for each component 
with those of HI. lBaan fc Irwinl jl995l) introduced a rotating 
disk model, in which the torus is confined in the roughly 
north-south orientation with an inclination about -20° w.r.t 
the CO(l-O) disk standing at P.A.= 15°. In this discussion, 



the jet axis connecting the nuclear radio sources of E and F 
makes an angle of about 45° with the confin ing torus. The 
sense of the rotation of the torus proposed in lBaan fc Irwinl 
( 1995) is consistent with that of the CO disk/core and 
the edge-on H2O maser disk proposed in |Trotte ret 
lll99Sl) and iKondratkd fcOOJl . ISawada-Satoh et alJ fcOOcf 
proposed a counter-rotating disk at P.A. = 30° in order to 
explain the spatially resolved HI absorption. However, the 
directions and position angles of the torus and the CO and 
H2O maser disks are completely inconsistent w ith those 
of the torus model in ISawada-Satoh et alJ <l2000l) and the 
evidence for this counter-rotation is weak. 

On the other hand, a counter-rotating structure re- 
versed to the kiloparsec-scale CO (1-0) disk has been ob- 
served in OH absorption in a type 2 Seyfert galaxy 
NGC5793 (Hagiwa ra et alJ 120001. Similarly, two nuclear 
disks with radii ~ 100 pc embedded in the outer kpc-scale 
gas disk have been clearly resolved in CO(2-l) in the merg- 
ing nuclei of the ultra-luminous infrared galaxy Arp 220 
(IC4553). The rotational sense between these two disks is 
reverse d due to the counter r otation of the two nuclei them- 
selves llSakamoto et alJll999l) . 



4-2.2 Interpretations of the double OH peaks 

The OH absorption in the EVN spectrum shows double 
pe aks with a separation of 101.6 km s^ 1 in velocity (Table 
2). lBaan fc Irwinl dl995l) considered the association of OH1 
and OH2 with the radio twin peaks E and F (Fig. 7) with an 
OH velocity gradient between the two components along the 
jet axis of approximately 50 cos <j> km s _1 pc -1 . The pro- 
jected separation of the twin radio peaks E and F is ~ 25 mas 
or 2 pc and the projected angle ((f)) lies between the radio 
axis E-F and the plane of the OH gas motion. The value of 
the gradient would be 10 times larger than that of the CO(l- 
0) kpc-scale disk and the CO(1-0) nuclear core rotating 
in the north-south direction, unless cos <j) is unrealistically 
small. Consequently, it is not plausible to correlate the OH 
velocity difference between E and F and the internal velocity 
gradients of the CO(1-0) disk/core, although the direction 
of the OH velocity gradients agrees with that of the larger- 
scale CO (1-0) disk. This velocity gradient would also be 
about 10 times larger than those seen in th e merging galaxy 
Arp 220 ijMundell. Ferruit fc Pedlarll200ll) and the Seyfert 
galaxies Mrk231 and Mrk27 3 ijKlockner. Baan. fc OarrettJ 
l2003t iKlockner fc Baanll2003l) . where HI absorption or OH 
maser emission reveal a rotating molecular torus with an in- 
ner radius of several tens of parsecs from the central engine. 
Considering the available evidence, the double peaks must 
arise from two kinematically independent systems in the nu- 
clear region. The direction of the velocity gradient of the sys- 
temic OH2 is in good agreement with that of the CO and 
H2O maser disks. In addition, the velocity range of OH2 
nearly coincides with the systemic velocity of NGC 3079. 
We know that the HI absorption at the systemic veloc- 
ity is seen against the wh ole nuclear continuum at 21-cm 
JSawada-Satoh et alJ I200H) . This can be accounted for by 
clumpy gas in a kiloparsec-scale disk. Accordingly, the sys- 
temic OH2 probes gaseous components in the CO disk, and 
particularly the inner disk on scales of 10-100 pc. 
On the other hand, the velocity field of the blue-shifted 
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QHl is very di f ferent from that of OH2 and of the CO disk. 
iBaan fe Irwinl jl 9951) argued that a foreground and an ex- 
panding shell driven by the nuclear superbubble could ex- 
plain the large blue-shift of the OHf centroid velocity, while 
a receding shell may account for the other weak red-shifted 
OH absorption component around Vlsr = 1260 km s _1 , as 
observed with VLA. In the starburst galaxy M82, distinct 
outflow components of molecular gas were discovered that 
extend over 500 pc above the plane of the disk and lie along 
the minor axis of the galaxy (e.g. Nakai et al. 1987). 
It is uncertain whether or not the outflows traced by OH1 
are associated with the (possible) starburst-related com- 
ponents E and F themselves because the velocity gradi- 
ent of OH1 does not align with this E-F axis (Fig. 7). It 
should be noted that the velocity range of OH1 (Vlsr = 
956-1050 km s _1 ) overlaps with that of several blue-shifted 
H2O features peaki ng at velocities Vls r = 1012, 1018, 1034, 
and 1035 km s" 1 ijTrotter et alJll998h . Although these fea- 
tures have a velocity signature similar to that of the OH1 
absorption component, there is no evidence for a physical as- 
sociation. Some of the H2O maser components outside the 
edge-on masering disk of the Circinus galaxy are interpreted 
as molecu lar outflow componen ts, ejected from the edge-on 
structure ijOreenhill et al.ll2003h . 

The two OH emission featu res identified at Vls r = 910 
km s -1 and 1230 km s _1 in lBaan fc Irwinl lll995l) and de- 
tected weakly in our data may be accounted for by the blue- 
shifted 1667 MHz and 1665 MHz OH features that lie at ve- 
locities just below the OH1 outflow components. They may 
arise in the shocked molecular components foreground to the 
outflow structure. In the above picture, an association with 
molecular outflows can explain the blue-shifted OH1, the 
weak OH emission, and also possibly the H2O maser com- 
ponents. However, this last association cannot be quanti- 
fied due to difficulties in comparing the distribution of weak 
H2O masers and the OH absorption on very different angu- 
lar scales. 



5 CONCLUSIONS AND SUMMARY 

The broad OH absorption towards the nuclear continuum 
source in NGC 3079 was imaged using VLBI techniques. The 
EVN observations reveal two kinematically independent OH 
absorption components, where OH1 is the blue-shifted ab- 
sorption, and OH2 is the absorption at the systemic velocity 
of the galaxy. The understanding of the kinematics of these 
components is limited by the fact that the 18-cm background 
continuum and the OH components have not been spatially 
resolved. 

1. The OH2 component shows a distinct velocity gradient 
in roughly north-to-south direction viewed against the un- 
resolved nuclear continuum at 18-cm consisting of compo- 
nents E and F. This velocity gradient of the systemic OH2 
is consistent with that of the clumpy gas component in the 
kpc-scale CO disk (P.A.= 15°). 

2. The blue-shifted component OH1 appears to be asso- 
ciated with molecular outflows and is possibly associated 
with the nuclear super-bubble due to circum-nuclear star- 
burst activity at components E and F and possibly at other 
locations. The velocity gradient of the OH1 component is 
not yet understood butis almost reversed from that of the 



systemic component OH2. The weak OH emission features 
marginally detected in our EVN observation may also be 
associated with shocked gas in these nuclear outflows. 
3. The H2O masers in the galaxy have been attributed to 
a compact disk structure around the true nucleus located 
west of the radio continuum component C. The OH1 veloc- 
ity range is similar to that of the blue-shifted H2O maser 
components and there may be some relation. 

The combined knowledge of our OH data, of previous 
VLBI data on the radio continuum, of H2O masers, and of 
the HI absorption on scales of 100 pc down to sub-parsec 
scales, provides a clear and consistent interpretation of the 
various spectral components seen at the nucleus. VLBI ob- 
servations at even higher angular resolution and high sensi- 
tivity could resolve the nature of the OH and HI components 
in relation to H2O maser components in the nuclear region 
of NGC 3079. 
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Figure 6. Comparison of the intensity distribution of the OH 
double peaks OH1 and OH2, integrated over 16 velocity channels 
(corresponding to ~ 180 km s _1 ) for each peak. Contour levels 
run from -2 to -1.2 mjy per beam in steps of 0.2 mjy per beam. 
The synthesized beam is included at the left corner of each plot. 
left: 1667 MHz, right: 1665 MHz. The dashed and solid lines show 
OH1 (blue-shifted absorption) and OH2 (near systemic velocity), 
respectively. 
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Figure 7. Schematic views of the nuclear region of NGC 3079. The sizes of components are drawn roughly in scale. Oval indicates 
the outline of the unresolved continuum nucleus in Fig. 2. Upward thin arrow denotes the direction of OH1 velocity gradient against E 
and F, both of them are not resolved in our EVN observations. Downward thick arrow shows the sense of the OH2 velocity gradient 
that is roughly consistent with the kpc-scale CO rotation, indicated by thin dashed arrows. A thin dashed lin e den otes the axis of 
the kiloparsec-scale super-bubble consisting of wide-angle outflows ejected from the nuclear region ICecil et al.ll200ll) . Inset displays 
the innermost nuclear region of NGC 3079, where three radio sources A, B, and C are aligned a long P.A.= 55°. The weak component 
C may be shaded by the inclined torus or disk traced by the H2O masers jTrot ter et alJll99aV Components E and F are visible at 
lower frequencies at from 1.7 to 5.0 GHz (Kondratko 2003; Middclbej^_et_alJ|2003j) and may result from circum - nuclea r starburst. The 
component labels are adopted from the nomenclature defined in llrwin fc SeaquisJ ll98°fl and iMiddelberg et ail 120031) . Components E 
and A show similar sp ectral characteristics, t hat are different from other components. The spectra of E and A are inverted at 2.3 and 
5.0 GHz, respectively iMiddelberg et al.ll2003l) . 



Table 2. Gaussian-fitted parameters of the OH absorption and emission spectra. VLA data are from Baan & Irwin (1995). Vc refers a 
line-peak velocity. Sp and AV are peak flux density and linewidth (FWHM), respectively. The EVN measurements of the OH emission 
arc tentative. 
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Table 3. Parameters of the 1667 MHz OH absorption. Global- VLBI data are from Sawada-Satoh et al. (2001) and VLA data are from 
Baan & Irwin (1995). 



EVN Global- VLBI VLA (A Configuration) 

OH1 (Blue-shift) OH2 (Systemic) OH1 OH2 OH1 OH2 



Optical depth (r) 0.46 ± 0.03 0.31 ± 0.02 0.73 ± 0.17 0.79 ± 0.13 1.03 1.28 

Nqh/Tex (K - 1 cm" 2 ) 4.3 x 10 15 5.3 x 10 15 9.2 x 10 1B 1.8 x 10 17 2.2 x 10 16 2.5 x 10 



